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We study the Mott insulator-superfluid transition of ultracold bosonic atoms in a two-dimensional 
square optical lattice in the presence of a synthetic magnetic field with p/q (p and q being co-prime 
integers) flux quanta passing through each lattice plaquette. We show that on approach to the 
transition from the Mott side, the momentum distribution of the bosons exhibits q precursor peaks 
within the first magnetic Brillouin zone. We also provide an effective theory for the transition and 
show that it involves q interacting boson fields. We construct, from a mean- field analysis of this 
effective theory, the superfluid ground states near the transition and compute, for g = 2, 3, both the 
gapped and the gap less collective modes of these states. We suggest experiments to test our theory. 

PACS numbers: 74.45+c, 74.78.Na 
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The physics of ultracold bosonic atoms in an opti- 
cal lattice can be well described by the Bose-Hubbard 
model [HQ. In fact, experiments on the Mott insulators- 
superfluid (MI-SF) transitions of such bosonic atoms 
in two-dimensional (2D) optical lattices Q is found to 
agree with predictions of theoretical studies of the Bose- 
Hubbard model quite accurately [4-6]. More recently, 
several experiments have successfully generated time- or 
space- dependent effective vector potentials for these neu- 
tral bosonic atoms by creating temporally or spatially 
dependent optical coupling between their internal states 
[1, [s[ . Such a generation of synthetic space-dependent 
vector potential and hence magnetic fields is complemen- 
tary to the conventional rotation technique [lO] . Several 
theoretical studies have also been carried on the prop- 
erties of the bosons in an optical lattice in the presence 
of an effective magnetic field fllj . In particular, the MI- 
SF phase boundary has been computed using mean-field 
theory [12| and excitation energy calculation u sing a per- 
tubative expansion in the hopping parameter [13|- How- 
ever, experimentally relevant issues such as the momen- 
tum distribution of the bosons in the Mott phase, the 
critical theory of the MI-SF transition, and the nature 
of the superfluid ground states and collective modes near 
criticality have not been addressed so far. 

In this letter, we present a theory of the MI-SF tran- 
sition for ultracold bosons in a 2D square optical lattice 
with commensurate filling no and in the presence of a syn- 
thetic vector potential corresponding to p/q (p and q are 
co-prime integers) flux quanta per plaquette of the lattice 
which addresses all of the above-mentioned issues. The 
novel results of our work which have not been addressed 
in earlier studies are as follows. First, using a strong- 
coupling RPA theory for the bosons [5] , we provide an an- 
alytical formula for their momentum distribution in the 
Mott phase and show that it develops q precursor peaks 
on approach to the MI-SF transition. Second, based on 
both the microscopic strong-coupling theory and a sym- 



metry analysis, we construct the critical field theory for 
the transition and show that it necessarily involves q cou- 
pled boson fields [14]. Third, using a mean- field analy- 
sis of this effective theory, we find the superfluid ground 
state to which the transition takes place and chart out the 
corresponding spatial patterns of the superfluid density. 
Finally, we compute the collective modes of the superfluid 
phase for q = 2,3, explicitly demonstrating the nature of 
both the gapped and gapless collective modes near the 
transition, and provide analytical expressions for their 
masses and group velocities in terms of microscopic pa- 
rameters of the theory. We suggest realistic experiments 
which can verify specific predictions of our theory. 

The Hamiltonian of a system of bosons in the presence 
of an optical lattice and a synthetic magnetic field is given 
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where J rr / = —Jex.p(—iq*f r A* • dl/hc), if r and r' 
are nearest neighboring sites and zero otherwise, A* = 
B*(Q,x) is the synthetic vector potential, q*(B*) is the 
effective charge (magnetic field) for the bosons, J is the 
hopping amplitude determined by the depth of the op- 
tical lattice, and the value of q*B* can be controlled by 
varying the detuning between the hyperfine states of the 
bosonic atoms |9j. Here /i is the chemical potential, U 
is the on-site Hubbard interaction, and b r (n r = blb r ) is 
the boson annihilation (density) operator. In the rest of 
this work, we consider the magnetic field to correspond to 
p/q flux quanta through the lattice: q*B*a 2 /he — 2itp/q, 
and set the lattice spacing a, h, and c to unity. 

The effect of the magnetic field manifests itself in the 
first term of Eq. [T] and thus vanishes in the local limit 
(J = 0). In this limit the boson Green function at 
T = can be exactly computed Q. 0. 0]: Go{iu n ) = 
(no + l)(iWn — Ep) -1 — no(iu n + Eh) -1 . Here uo n denote 
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FIG. 1: Color online) The MI-SF phase boundary for q = 2. 
The red dot indicates the value of J and p for which n(k) in 
the left panel Fig. [2] has been plotted. 



bosonic Matsubara frequencies and Eh = \i — U(no — 
l)(E p = — fi + Uno) are the energy cost of adding a hole 
(particle) to the Mott state. To address the effects of the 
hopping term, we write down the coherent state path 
integral corresponding to T~L: Z = J Di/jDip* exp(— S) 

where S = tfdT[fc r fc(T)d T $ r (T)+H$\4])], t is 
the imaginary time, j3 = l/fc^T is the inverse tem- 
perature (T), and ks is the Boltzman constant. Fol- 
lowing Ref. |5j, we then decouple the hopping term by 
two successive Hubbard- Stratonovitch transformations, 
integrate out the original boson and the first Hubbard- 
Stratonovitch fields, and obtain the final form of the 
strong-coupling effective action S e fi = So + S\ 

S = ^*(io; n ,k)[-Go 1 (iu n )I + J q (k)]ip q (icj n ,k), 



(2) 



Si = 9l*j*dr j ^ 2 r|^*(r,r)^(r,r)| 2 



where ijj q denotes the ^-component auxiliary field intro- 
duced through the second Hubbard- Stratonovich trans- 
formation and have the same correlation functions as the 
original boson fields $ [5], / k = (1//3) £ Wti / d 2 k/(2ir) 2 , 
I denotes the unit matrix, and g > is the static limit 
of the exact two-particle vertex function of the bosons in 
the local limit @. Here J,(k) is a q x q dimensional tridi- 
agonal hermitian matrix whose upper off-diagonal [diag- 
onal] elements are — Jex.p(—ik y ) [—2Jcos(k x +2ira/q)], 
with a = 0, 1, ..g — 1. It is well-known that J q (k) has q 
eigenvalues e q (k) within the first magnetic Brillouin zone 
(— 7r < k y < 7r, — 7r/q < k x < n/q) which are g-fold de- 
generate. In particular, the lowest eigenvalue e q n (k) has 
q degenerate minima at Q a = (0,27ra/g) [l5| . Note that 
So reproduces correct bosons propagator both in the lo- 
cal (J = 0) and the non-interacting (U = 0) limits. Also, 
since Gq 1 is independent of momenta, finding the bo- 
son Green function G(zcj n ,k) = [— Gq 1 (iuj n )I + J g (k)] _1 
amounts to inverting J g (k). 





FIG. 2: (Color online) Plot of n(k) for q = 2 (left panel) 
and q = 4 ( right panel) at /x/C7 = 0.414 and J/J c = 0.95 
indicating the precursor peaks in the Mott phase. 



The critical hopping J c for the MI-SF transition as a 
function of \i can be determined from the condition Q 

r q = -Go 1 (iw„ = 0) + e™(k = Q a )=0. (3) 

The MI-SF phase boundary so obtained is shown in Fig.[T] 
for q = 2 and agrees qualitatively with those obtained us- 
ing mean-field theory [12] and J/U expansion [13[. Note 
that J c remains same for all Q a due to the g-fold degen- 
eracy of e™(k). 

The consequence of the q fold degeneracy of e q becomes 
evident in the momentum distribution of the bosons in 
the Mott phase, which, at T = 0, is given by n(k) = 
limT-K)(l//3) ^Zojn TrG(zo; n , k). After some straightfor- 
ward algebra, one obtains 



i(k) 
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where S/i — ji — U(uq — 1/2), p = (no + 1/2) and 
£«±(k) = -Sfi + e»(k)/2 ± x /e«(k)2 + 4e«(k)Up + m/2 



denotes the position of the poles of G(k, iuj n ) in the Mott 
phase. Note that E q can also be obtained from a time- 
dependent variational method [l6| . 

Eq. [H is a central result of this work and generalizes its 
counterpart in Ref. 5 in the presence of a magnetic field. 
The peaks of n(k) occur when E q + (k) — E q ~ (k) becomes 
small near the MI-SF transition. The degeneracy of e"(k) 
and hence ^^(k) ensures that this happens at q points 
in the first Brillouin zone leading to q precursor peaks 
in n(k) at k = Q a . This is demonstrated in Fig. [2] 
for q = 2 and q = 4. Note that the positions of these 
peaks in the Brillouin zone depend on the specific form 
of the vector potential realized in the experiments; for 
symmetric vector potentials generated by rotation they 
would appear at (na/q^Tra/q). However, their number 
depends only on p/q and the lattice geometry. 

At J c , the MI-SF transition occurs since the energy gap 
to addition of particles and/or holes to the Mott state 
vanishes. In contrast to standard superfluid-insulator 
transition 0-0, the presence of q degenerate minima at 
k = Q a necessitates the corresponding Landau-Ginzburg 
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FIG. 3: (Color online) Plot of the normalized superfluid den- 
sity p^ for q — 2 (left panel) and pi 3 ^ for q = 3 (right panel). 

theory to be constructed out of q low-energy fluctuating 
fields (f) a (r,t) around these minima: 

^(r,i) = Yl X?WM), (5) 

a=0..q-l 

where Xq( T ) denotes the eigenvectors of J q (Q a ) in real 
space, and we have Wick-rotated to real time. The 
quadratic part of the Landau-Ginzburg theory, obtained 
by expanding So (Eq. [2]) about the minima, is given by 

So = fd 2 rdt ^(r,r)[K d 2 ^iKid t 

^ a=0..q-l 

+r,-^(^ + ^)]^(r,r), (6) 

where K = l/2d 2 G^ /dou%= = n (n + l)U 2 /(fi + 
Z7) 3 , Kx = 5G?o V^U=o = 1 - Mno + 1)^70" + U) 2 , 
and v 2 = V£e™(k)/2 with v\ = J/y/2. At the tip of 

the Mott lobe, where /i = /i t ip = ^(ynoK + l) — 1), 
K\ = 0. Thus we have a critical theory with dynamical 
critical exponent z = 1. Away from the tip, i^i ^ 
rendering z = 2 [2|. 

The most general quartic Landau-Ginzburg action 
in terms of q bosonic fields which is allowed by in- 
variance under projective symmetry group (PSG) of 
the underlying square lattice has been obtained in 
Ref. [14]. The elements of the PSG for the square 
lattice include translation along x and rotation 
by 7r/2 about z axis, and reflections about x and 
y axes. The transformation properties of (j) a fields 
under these operations are tabulated in Ref. 
The invariant quartic action so obtained is given 
by Si = / ' dVdt X^ )7=0 ^ 

where Y^ = Y~ a ~^ = Y^~^ = T^~ 2 ^-^ and sums 
over integers a, /?, and 7 are taken modulo q. Eq. [6] along 
with Si has been analyzed in details in Ref. [14(. How- 
ever, the lack of microscopic knowledge of Y^ did not 
allow identification of the exact ground state of <Si; only 
possible symmetry- allowed ground states were charted. 

Here, taking advantage of the microscopic knowledge 
of g and Xr ^ we determine the exact superfluid state to 
which the transition takes place. This is done by sub- 
stituting of Eq. [5] in Eq. [2] followed by coarse-graining of 



the resultant action which involves replacing Xr by its 
sum over a lattice sites: / d 2 rdtL 1 [Xq(r)]L 2 [<p a {r,t)] -> 
{W)El:* Li[x%(r)]}Jd 2 rdtL 2 [r(r,t)} = c J d 2 r 
dt L2[</> a (r, t)]. Here Li and L 2 denotes arbitrary fourth 
order polynomial functions and the coarse-graining pro- 
cedure is applicable due to the natural separation of 
scale between the spatial variation of Xq ( r ) an d (/> a (r, t). 
The effective action so obtained is then compared to S\ 
to obtain Y^ . Finally, we minimize the resultant ac- 
tion at the mean-field level and obtain the superfluid 
ground state near the MI-SF transition. This procedure 
is most easily demonstrated for q = 2. Here, ^(k) = 
— 2Jy/cos 2 (k x ) + cos 2 (k y ) leading to two minima at 
(k x ,k y ) = (0, 0) and (0, 7r) with eigenfunctions X2( r ) = 
(l + v / 2+exp(z7rx))/v / 4 + 2y/2 and x 2 ( r ) = exp(z7n/)(l + 
\/2 — exp(i7rx))/ y4 + 2\/2. Putting these values in 
Eq. [5] and Eq. [2j the coarse-grained effective action 
reads S q = 2 = 1/8 / d 2 rdt[3g(\^° (r, t)\ 2 + ^(r, t)\ 2 ) 2 + 
g(0 o *(r,t)^(r,t) - ^(r, t)0°(r, t)) 2 ]. Comparing S q = 2 
with Si for q = 2, we find r°° = 3g/2 and T^ = g/2. 
A mean-field analysis then yields the superfluid ground 
state: (0°(r,t)) = 0° = icj) 1 = (^(r,*)). The renor- 
malized superfluid density can be obtained by using 
p { s\v) = |^ F (r)| 2 /|^ F (0)| 2 where ^f ¥ is obtained by 
substituting ((/> 0,1 (r, t)) in Eq. [51 Analogous procedure 
carried out for q = 3 yields the superfluid ground state: 
(0°(r,t)) ^ 0, ((/> a * (r,t)) = 0. The resultant plots of 
Ps^(r), shown in Fig. |4]for q = 2, and q = 3, display 2 
and 3 sublattice patterns respectively. We note that the 
procedure mentioned above constitutes a general method 
for obtaining the superfluid ground state and density near 
the MI-SF critical point for any q. 

Finally, we compute the collective modes of the su- 
perfluid ground state near the transition. First we 

(2) 

consider the case q — 2 and rewrite S^ J in terms of 
a linear combination of the <j) a fields: £°M = (0° -f 
[— ^cj) 1 )/^. The quartic action becomes S e ^~ 2 = 
1/8 J d 2 rdt[U\e(r,t)\ 2 + I^M)! 2 ) 2 - 5(|£°M)| 2 - 
^(r, t)| 2 ) 2 ] so that the superfluid ground state corre- 
sponds to condensation of £°. The quadratic action can 
be written as S = J d 2 rdt^ q ~J ^(r,t)[-GQ 1 (uj) - 
c 2 + ^|k| 2 ]e*(r,£), with c 2 = -e 2 (k = 0) = 2^/2 J . Us- 
ing these actions, and carrying out a straightforward lin- 
earization g°( r,t) = f + 5l°(r,t) and^(r,t) = S^(r,t), 
where £° = y/2\r2\/g, we find that there are four collec- 
tive modes. Two of these correspond to the condensed 
field £° and £ *, and have dispersions 

4 1} = (±B 2 (k)/2+ [Bl(k)/4-C 2 (k)] 1/2 ) 1 /2, ( 7 ) 

where B 2 (k) = [2<5/x - A 2 (k)} 2 + 2a (w||k| 2 - r 2 ) - r 2 2 , 
C 2 (k) = alv 2 2 \k\ 2 {v 2 2 \k\ 2 - 2r 2 ), a = (U + fi), and 
A 2 (k) = —c 2 + w 2 |k| 2 — 2r 2 . At low wave-vector, lo^ 
is gapped with a mass mi 

ear dispersion with velocity vq = ^2<^o\/^l r 2|/ m i- The 
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FIG. 4: (Color online) Top left, right, and bottom left panels: 
Plot of mi, 2, 3 as a function of \i for J/J c = 1 (solid black 
line), 1.2 (red dotted line) and 1.5 (blue dashed line). Bottom 
right panel: Plot of vg/v2 vs \i for J/J c = 1.01 (black solid 
line), 1.2 (red dotted line) and 1.5 (blue dashed line). U = 1, 
no = 1, <? = 2, and /xtip = 0.414 for all the plots. 



For experimental verification of our theory, we suggest 
measurement of n(k) for the bosons in the Mott phase 
near the transition as done earlier in Ref. Q for 2D opti- 
cal lattices without the synthetic magnetic field. This 
distribution is predicted to display q precursor peaks. 
The collective modes in the superfluid phase can also 
be directly probed and compared to the theory by stan- 
dard lattice modulation experiments [17j and response 
functions measurement by Bragg spectroscopy [Hj . 

In conclusion, we have analyzed the MI-SF transition 
of ultracold bosons in a 2D optical lattice in the presence 
of a synthetic magnetic field. We have demonstrated the 
presence of q precursor peaks in their momentum distri- 
bution near the MI-SF transition, provided a critical field 
theory for the transition, analyzed this theory to predict 
the ground state and the collective modes of the bosons 
in the superfluid phase, and suggested experiments to 
probe our theory. K.S. thanks DST, India for financial 
support under Project No. SR/S2/CMP-001/2009. 



other two modes, which correspond to the non-condensed 
field x 1 an d X 1 *? nave dispersions 

4 2) = ±D 2 (k)/2 + [D 2 (k) 2 /4 + a (\r 2 \/2 + v 2 \k\ 2 ))^ 2 , 

where D 2 {k) = -(25fi + c 2 - ^|k| 2 + r 2 /2). Both 
these modes are gapped in the superfluid phase with 
masses m 2[3] = +[-] J D 2 (0)/2 + y/D 2 (0) 2 /4~+ a \r 2 \/2. 
The masses rai ? 2,3 and the velocity vq of these modes, 
plotted as a function \i in Fig. |4] for several representa- 
tive values of J/J c , displays the following characteristics. 
At \i = /itip and J = J c , where 25 \i — —c 2 rendering 
B 2 (0) = and D 2 (0) = 0, all the modes become gapless 
with uj ~ |k| dispersion. Also at fi ^ /i t i P , one of the two 
modes always remain gapless at J = J c with uo ~ |k| 2 
dispersion. The velocity vq at J = J c , is non-zero only 
at fi = /itip; thus it shows a peak at /i t i P for J close to J c . 
We emphasize that our theory specifies vq and mi ( 2,3 in 
terms of the parameters of the Bose-Hubbard model. 

For q = 3 only condense, and the corre- 
sponding collective modes are given by Eq. with 
c 2 ,v 2 ,r 2 c 3 ,v 3 ,r 3 (where c 3 = -e^(0)). This 
leads to similar gapped and a gapless mode with lin- 
ear dispersion as for q = 2. However, the disper- 
sion of the non-condensed modes are different. The 
effective action S e ^~ turns out to be 0(3) sym- 

metric: <|= 3 ~ / d 2 rdt(E a=0 .. 2 l£ a (r, t)\ 2 ) 2 lead- 
ing to two doubly- degenerate non-condensed modes 

4 3) = (±233 (k) + ^/D 3 (k) 2 + v 3 \k\ 2 )/2, where D 3 (k) = 
— (25{i + C3 — v 2 \k.\ 2 ). Thus there are two gapped and two 
gapless modes with u ~ |k| 2 . These two modes become 
gapless due to the 0(3) symmetric form of S^ -3 . For 
q > 3, there are in general 2q collective modes, and we 
have left their analysis as a subject of future study. 
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